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PARABOLIC  EQUATION  STARTING  FIELD 
FOR  A  PROLATE  SPHEROID  SOURCE 


INTRODUCTION 

Parabolic  equations  are  used  to  calculate  the  propagation  of  acoustic 
fields  over  long  distances  in  waveguides.  To  use  parabolic  equations  one 
must  specify  initial  values.  In  the  case  of  sound  propagation  in  the  ocean, 
the  inital  values  required  are  the  so-called  "starting  field"  which  must  be 
given  from  the  surface  to  the  bottom  of  the  ocean  completely  surounding  the 
source  of  sound.  Once  this  is  given,  a  numerical  parabolic  equation  model 
can  be  used  to  advance  this  field  one  horizontal  step  at  a  time  to  any 
desired  range.  We  are  interested  in  calculating  such  fields  at  rather  close 
ranges.  These  short  ranges  imply  that  the  source  should  not  be  considered 
to  be  a  point  source.  For  this  paper  we  model  very  special  sources  as 
vibrating  prolate  spheroids,  such  sources  being  cigar  shaped.  We  examine 
this  special  problem  to  give  an  idea  of  the  effects  on  the  starting  field, 
and  hence  on  the  acoustic  field  when  the  source  is  not  a  point  source. 


SEPARATION  IN  PROLATE  SPHEROIDAL  COORDINATES 

The  sound  field  in  a  uniform  medium  is  governed  by  the  wave  equation 
given  by 

V2*  -  -2  0tt.  (1) 

c 

We  want  to  solve  this  equation  in  a  uniform  halfspace  that  contains  a 
vibrating  prolate  spheroid.  We  follow  the  usual  approach  of  eigenvalue 
expansion  [ref.  1-3].  Instead  of  the  normal  Euclidean  coordinate  system 
(x,  y,  z)  a  prolate  spheroidal  coordinate  system  (^,  8,  4>'  is  used  where  the 
two  systems  are  related  by 

x  -  a/2  sinh(fs)  sin(8 )cos  (<j>) ,  (2) 

y  -  a/2  sinh(p)sin($)sin(^)  and  (3) 

z  -  a/2  cosh(M)cos(ff ) .  (4) 

If  we  write  cosh(^)  -  ?  and  cos($)  -  rj ,  then  we  get 


-fJ^  f 


n  cos(^) 


(5) 
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y 


sin(^) 


and 


(6) 


where  the  ranges  of  i|  and  <  are  given  by  latfs*.  -lSn^l  and  0s^<2*. 
Prolate  spheroids  are  defined  by  n  -  constant  or  {  -  constant;  whereas  B- 
constant  or  *?  -  constant  give  hyperboloids  of  two  sheets  and  4>  -  constant 
gives  the  half  planes  through  the  z*axis  (See  Figure  1). 


Figure  1 

PROLATE  SPHEROIDAL  RADIATOR 


If  we  substitute 


4>  -  S(-,)R(f)JWe'iut 

into  equation  (1)  we  get  the  ordinary  differential  equations 


where  A  and  m  are  separation  constants,  and  h-ka/2  where  k-u>/c . 
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As  h-*0,  the  spheroidal  coordinates  go  to  spherical  co-ordinates  and  the 
S-functions  go  to  associated  Legendre  polynomials  P^(?j).  When  h*0  we 
denote  the  dependence  of  S-functions  on  i ,  m  and  h  and  write  them  as 


S(tj) 

or 

S(i») 


S(h,t),ro,i)  “  £  d(h,m,  i  ,  2n)P™+2n(»?) ,  when  m-i  is  0  or  even 
n-0  *“ 

00 

S(h,fj,m,i)  -  X  d(h,m,i  •  2n+1)p^+2n+l when  is  odd  • 
n-0 


(12) 


(1?) 


The  standard  notation  in  the  literature  [refs.  1,2,4]  is  given  as 
S(tj)  -  E'  d(h,m,  i,n)P^+n(»j)  where  the  '  denotes  the  proper  choice  of  (12)  or 
(13)  depending  on  the  parity  of  (m-i). 


The  radial  solutions  R(£)  show  a 
functions  S (rj) .  The  two  solutions  ' 
of  spherical  Bessel  functions  as  follows 


structure  similar  to  the  angular 
(2)  ° 
and  Rv  '  can  be  expressed  in  terms 


r  2 

-»m/2 

00 

R(1ta, inl¬ 

r- 

_1 

}_  a(h.m.i.n)jm+n(hO 

*2 

k, 

n-0 

and 

9 

^m/2 

oo 

_ _ _  • 

R(2ta,m,i)- 

r- 

_1 

a(h,m,i,n)nit+n(hO 

n-0 

C2 

4 

(14) 


(15) 


For  expressing  complex  radial  solutions  propagating  outward  from  the  prolate 
spheroid  in  terms  of  spherical  Hankel  functions  we  use 

r(3)  _  R(l)  +  iR( 2 )  (16) 

because  it  satisfies  the  Sommerfeld  radiation  condition. 


A  VIBRATING  PROLATE  SPHEROIDAL  SOURCE 

We  consider  a  prolate  spheroid  whose  surface  vibrating  in  its 

fundamental  mode  with  a  velocity 

v  -  S  (h ,  ij  ,  0 , 0)  (£q  -  rj2r1/2  e*iwt.  (17) 

We  seek  the  corresponding  pressure  field  in  the  form 

p  -  A  S(h.n.O.O)  R(h?,0,0)  e*iwt  ,  (18) 

Where  A  is  a  constant.  The  balance  of  force,  mass  and  accel ’ration  is  [ref. 

6] 
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„  3v 

yP  ~  ~p°  Jt  ’ 

which  expanded  in  vector  form  is  [ref.  6] 

1  1 


(19) 


82 


c2  2  an  ’  [? 
s  o  *  n  s  o 


=  iE  o 
2  ac  ’ 


iwpt 


°.  ri - —  S(h.9fOfO)e*lwt.O 


2 

,  ^0 


(20) 


The  first  two  components  evaluated  at  the  surface  of  the  prolate  spheroid 
give 


2  2 

£0  ■  n 


A  S'(h,»7,0,0)R(h£0l0,0)e 


-  iwt 


-  0 


(21) 


and 


hA 


2  2 
,  ?0  *  n 


S(h,i7,0,0)R'(h^0,0, 0)eiwC 


|  2  2 

i«o  -  n 


S(h,q,0,0)eiwt  (22) 


Equation  (21)  states  that  £0  must  be  a  zero  of  R(h£,0,0),  that  is  to  say,  £0 
must  be  an  eigenvalue.  Equation  (22)  is  the  equation  from  which  we  derive 
the  value  of  A 


A  - 


-  law 


■P  o 


2h  R(3)(h|o,0,0) 


(23) 


The  pressure  p  is  therefore  expressed  as 

iawpo  R(3)(h£,0,0)  „  „  „  -iwt 

P  “  — jr-  -  S(h,q  ,0,0)  e 

RU;(hCo,0,0) 

where  £0  is  some  eigenvalue  of  the  vibration  of  the  prolate  spheroid. 


(24) 


TRANSFORMATION  TO  THE  USUAL  OCEAN  COORDINATES 

The  prolate  spheroid  with  its  axis  of  symetry  on  the  z-axis  of  a 
rectangular  (x,y,z)  coordinate  system  is  transformed  to  a  horizontal  prolate 
spheroid  in  a  rectangular  (p.x.f)  coordinate  system  by  the  transformation 
(See  Figure  2) 


P  -  z  , 

X  m  y  and 

r  -  x  +  rc 


(25) 

(26) 
(27) 


The  transformation  takes  the  center  of  the  prolate  spheroid  in  the  (x,y,z) 
coordinate  system  to  the  point  (0,0, f0)  in  the  (p.x.D  coordinate  system  and 
the  axis  of  symmetry  is  now  parallel  to  the  p  axis. 


u 


Figure  2 


TRANSFORMATION  OF  COORDINATES 
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If  we  solve  the  equations  (2) -(A)  for  £  and  »?  we  get 


£  -  -  x2+y2+(z+a/2)2  +  -  |x2+y2+ 

a 


+y  +(z-a/2) ' 


i,  -  ^  x2+y2+(z+a/2)2  -  *  x2+y2+(z-a/2) 2  . 


When  x,  y  and  z  are  expressed  in  terms  of  p,  x  an&  f  we  obtain 
£  -  \  (f -f0)2+x2+(p+a/2)2  +  -  ,  (r-ro>2+x2+(p-a/2)2 


n  ~  ~  (r-r0>2+x2+(p+a/2)2  *  a  i  •fo)  +x?+(P"a/2) 


METHOD  OF  IMAGES 


To  generate  our  starting  field  we  use  the  method  of  images  involving 
two  identical  prolate  spheroidal  sources,  opposite  in  phase;  one  at  a  depth 
of  f0  and  another  at  -f0  (See  Figure  3).  This  allows  us  to  match  a  pressure 
release  condition  (p-0)  at  the  surface  of  the  ocean. 
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Figure  3 


The  tocal  pressure  is  the  superposition  of  the  field  from  each  individual 
source  which  is 

p  -  const  eiwt[R(3)(h£,0,0)S<h,f?.0,0)  -  R(3)  (h£  '  ,  0 , 0)S(h ,  r, '  ,  0 , 0)  ]  (32) 

where 

r-  2  J  (S+t0)7+X2  +  (p+a/2)7  + 

»»'-  \  J  (f+f0)2+^2+(p+a/2)3  * 

Finally,  we  obtain  the  starting  field  that  corresponds  to  our  special 
source,  a  prolate  spheroid  vibrating  in  its  fundamental  mode,  by  sampling 
the  expression  in  (32)  over  the  entire  depth  of  the  ocean  at  a  fixed  range. 
This  starting  field  will  vary  with  the  direction  from  the  source  because  the 
source  is  not  rotationally  symmetric  about  the  C  axis. 


2  ,  (c+ro)2+*2+(p'a/2)2  and 


(33) 


2  <r+$-0)2+^2+^-a/2): 


(34) 


6 


TD  8035 


REFERENCES 

1.  J.  A.  Stratton  et.  al . ,  Spheroidal  Wave  Functions,  Technology  Press  of 
M.I.T.  and  John  Wiley  &  Sons,  New  York,  1956. 

2.  C.  Flammer,  Spheroidal  Wave  Functions,  Stanford  University  Press, 

Stanford  CA,  1957. 

3.  S.  Hanish,  A  Treatise  on  Acoustic  Radiation,  Naval  Research  Laboratory, 
Washington,  D.C. ,  1981. 

4.  P.  M.  Morse  and  H.  Feshbach,  Methods  of  Mathematical  Physics,  McGraw- 
Hill,  New  York,  1953. 

5.  P.  M.  Morse  and  U.  Ingard,  Theoretical  Acoustics.  McGraw-Hill,  New 

York,  1968. 

6.  P.  M.  Moon,  and  D.  E.  Spencer,  Field  Theory  Handbook:  Including 

Coordinate  Systems,  Differential  Equations  and  Their  Solutions, 
Springer  Verlag,  New  York,  1971. 


7/S 

Reverse  Blank 


INITIAL  DISTRIBUTION  LIST 


Addressee  No.  of  Copies 

SACLANT  ASV  Research  Center  (Dr.  R.  L.  Martin,  Dr.  F.  c.  Jensen)  2 

NRL  (Dr.  Ralph  Baer,  Code  5160;  Orest  Dlachok;  Dr.  John  s. 

Perkins,  Code  5160;  Dr.  V.  A.  Kuperman;  Dr.  Michael  Porter)  5 

NOSC  (Code  8302;  Dr.  Hooter  Bucker,  Code  5311B;  David  Gordon; 

M.  A.  Pedersen)  4 

NORDA  (Dr.  Stanley  Chin-Bing;  Dr.  H.  F.  Werby;  R.  V.  McGirr  (1); 

G.  Morris  (1);  Edward  Chaika  (l);  Dr.  David  King)  6 

NCSC  (Dr.  Larry  Green,  Code  4120)  1 

CNO  (OP-096)  1 

NAVSEASYSCOM  (SEA-003,  -06R,  -63D,  -63R,  -63R-1,  -63R-13,  -53RA)  *7 

NAVAIRDEVCEN  1 

DTIC  2 

CNA  1 

Defense  Research  Establishment,  Pacific  (Dr.  David  Thomson, 

Dr.  Gary  Brooke,  Dr.  N.  Ross  Chapman)  3 


Defense  Research  Establishment,  Atlantic  (Dr.  Dale  D.  Ellis, 

Dr.  D.  Chapman) 

NOAA/ERL/WDL  (Michael  Jones) 

NQAA/AOML/OAL  (Dr.  David  Palmer,  Code  5122;  Dr.  John  Tsai) 

NOO,  NSTL  Station  (Nils  Paz,  Code  *7332;  Dr.  Lloyd  R.  Breslau, 

Code  9000;  Dr.  William  Jobst) 

ONR,  Arlington  (Robert  Obrochta,  Dr.  Richard  Lau, 

Dr.  Marshall  Orr,  Dr.  Raymond  Fitzgerald) 

ONR  (Codes  100,  200,  220,  222,  400,  422,  425AC,  480,  486) 

OCEANAV 

NAVOCEANO 

Weapons  System  Research  Laboratory  (Dr.  D.  Kewley) 

Georgia  Institute  of  Technology  (Prof.  William  P.  Ames, 

Prof.  Allan  D.  Pierce) 

University  of  Miami  (Prof.  Harry  DeFarrai,  Prof.  Frederick 
D.  Tapper t) 

Colorado  School  of  Mine  (Dr.  John  A.  DeSanto;  Prof.  L.  Fishman 

(Center  for  Wave  Phenomena,  Dept  of  Math;  Golden,  CO  80401) 
Science  Applications,  Inc.  (Dr.  Lewis  Dozier,  Dr.  John  Hanna, 

Dr.  C.  W.  Spofford) 

University  of  Toronto,  Canada  (Prof.  Kenneth  R.  Jackson) 

ODSI  Defense  System  Inc.  (Paul  Etter) 

APL/Unlversity  of  Texas  (Prof.  Robert  A.  Koch,  Susan  Payne) 
APL/Penn  State  (Dr.  S.  T.  McDaniel,  Prof.  Allan  D.  Pierce) 
Northwestern  University  (Prof.  G.  Kriegsmann) 

Columbia  University  (Prof.  John  T.  Kuo,  Dr.  Yu-Chiung  Teng) 


Admiralty  Underwater  Weapons  Establishment  (Dr.  Roy  G.  Levers)  1 
University  of  Rhode  Island  (Dr.  Larry  r.ayer,  Prof.  G.  R.  Verma)  2 
University  of  New  Orleans  (Prof.  J.  E.  Murphy)  1 
Daubin  Systems  Corp.  (Dr.  Lan  Nghiem-Phu)  1 
Bell  Telephone  Laboratory  (Richard  S.  Patton,  Denni  Seals)  2 
Polytechnic  in  Westchester  (Prof.  S.  Preiser)  1 
Yale  University  (Prof.  Martin  H.  Schultz)  4 


<N  CM  <N  00  t)*  O'  ■-<  ON  ^  (N  ON  ON  01  V  h  N  N  N 


INITIAL  DISTRIBUTION  LIST  (Cont’d) 


Rensselaer  Polytechnic  Institute  (Prof.  w.  L.  siegoann, 

Prof.  M.  J.  Jacobson)  6 
University  of  Massachusetts  (Prof.  Donald  P.  St.  Mary)  6 
Rutgers  University  (Prof.  Robert  Vlchnevetsky)  1 
Iowa  State  University  (Dr.  James  corones)  1 
FLOPETROL,  SCHLUMBERGER .  Prance  1 
Porschungsanstalt  der  Bundeswehr  fuer  Vasserschall  un  Qeophysik 

(Prof.  G.  Ziehm)  1 
University  of  Crete,  Greece  (Prof.  John  S.  Papadakis)  1 
Harvard  university  (Prof.  Allan  R.  Robinson,  Dr.  Wayne  Leslie, 

Dr.  Scott  Glenn,  LT  Jesse  Carmen)  4 
Massachusetts  Institute  of  Technology  (Prof.  A.  B.  Baggeroer)  1 
INO  (Dr.  Christopher  Mooers,  Dr.  David  Adamec,  Dr.  Michele 

Rienecker,  Robert  Willems,  Dr.  Chlng-Sang  Chul)  5 
The  Catholic  University  of  America  (Prof.  J.  J.  McCoy)  1 
The  Florida  State  University  (Prof.  E.  C.  Young,  Prof.  C.  Tam)  2 
Dr.  Michael  J.  Buckingham  (Royal  Aircraft  Establishment 

Farnborough  Hants  GU14  6TD  ENGLAND)  1 


